EPJ manuscript No. 

(will be inserted by the editor) 



o 
o 

> 
O 

o 

(N 



o 



> 
in 

O 

o 



o 
o 



X 



Breakdown of Luttinger's theorem in two-orbital Mott insulators 

A. Rosch 

Institute for Theoretical Physics, University of Cologne, 50937 Cologne, Germany 
Received: date / Revised version: date 

Abstract. An analysis of Luttinger's theorem shows that - contrary to recent claims - it is not valid for 
a generic Mott insulator. For a two-orbital Hubbard model with two electrons per site the crossover from 
a non-magnetic correlated insulating phase (Mott or Kondo insulator) to a band insulator is investigated. 
Mott insulating phases are characterized by poles of the self-energy and corresponding zeros in the Greens 
functions defining a "Luttinger surface" which is absent for band insulators. Nevertheless, the ground states 
of such insulators with two electrons per unit cell are adiabatically connected. 

PACS. 71.10.-w Theories and models of many-electron systems - 71.10.Fd Lattice fermion models (Hub- 
bard model, etc.) - 71.30.-|-h Metal-insulator transitions and other electronic transitions 



The basic quantity which defines a metal at low tem- 
peratures is the Fermi surface. Excitations of the Fermi 
surface are the basis of Fermi liquid theory. The concept 
of the Fermi surface also allows to distinguish different 
phases. Changes in the topology of the Fermi surfaces 
(e.g. vanishing bands) are therefore always associated with 
quantum phase transitions. 

Does an analog quantity exist in the case of a Mott in- 
sulator where strong interactions prohibit the formation 
of a Fermi surface? In a number of recent papers [UElElIll 
[5l[6] it has been argued that in such a situation the con- 
cept of a Fermi surface has to be replaced by the so-called 
"Luttinger surface". While the Fermi surface is defined 
by poles of the Greens function, the Luttinger surface is 
obtained from the zeros of the Greens function or, equiv- 
alently, from the poles of the self-energy. 

This identification is motivated by Luttinger's theo- 
rem [8], which states that at zero temperature, T = 0, 
the total density of electrons n can be obtained by sum- 
ming up all momenta where the real part of the Greens 
function evaluated at the Fermi energy is positive. For a 
multiband system the propagator is a matrix. Denoting 
the eigenvalues of this matrix by Gq,(w, p) and using that 
lmGa{0,p) — 0, Luttinger's theorem (also called Lut- 
tinger's sum rule) takes the form 



n = 2 



Ga(w = 0,p)> 



(2^)^ 



(1) 



It claims that the total density of electrons is fixed by the 
volume enclosed by a surface where the zero-frequency 
propagator G(0,p) changes sign (alternative versions of 
the theorem are briefly discussed in Sec. [4]). As emphasized 
recently by Dzyaloshinskii [IJ and by Essler and Tsve- 
lik [5], there are two ways how such a sign change can 
happen. At a Fermi surface, E{p) = pi, there is a pole. 



G(w, p) w ^-(g(p)-p) ' where Zp and Ep are the weight 
and dispersion of the quasi particle. Alternatively, the 
propagator can change its sign going through zero instead 
of infinity, G{uj = 0, p) = 0. The latter condition defines 
the "Luttinger surface" which is realized in Mott insula- 
tors [l,2j3j. This concept was for example used by Yang, 
Rice and Zhang [4] to build up a phenomenological the- 
ory of the pseudogap state. In a recent preprint, Stanescu, 
Phillips and Choy |5] tried to explore the role of Luttinger 
surfaces for (doped) Mott insulators. Recently, Ortloff, 
Balzer and Potthoff |[7j investigated under what conditions 
various (non-perturbative) approximation schemes lead to 
a violation of the Luttinger theorem. 



This motivates us to investigate the question whether 
the concept of such a Luttinger surface is as general and 
robust as the Fermi surface. Can one conclude that two 
states with a different topology of Luttinger surfaces have 
to be separated by a quantum phase transition? In Ref . [3] , 
Konik, Rice and Tsvelik have for example constructed a 
doped spin liquid (by coupling Mott-insulating Hubbard 
ladders) which is characterized both by a Luttinger surface 
and by the Fermi surfaces of small particle and hole pock- 
ets. Here the question arises whether this seemingly exotic 
state of matter is adiabatically connected to a weakly in- 
teracting system where Luttinger surfaces are absent. 



In the following we will first show that Luttinger's the- 
orem (HI) is not valid for a generic Mott insulator using an 
explicit counter example. We will then discuss which as- 
sumptions underlying its proof may not be fulfilled and 
finally investigate implications for the question whether 
various insulating states are adiabatically connected. 
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1 Model 

For definiteness, the following two-band Hubbard model 
is considered 



H - Hkin + ^ ffloc (2) 
i 

Hkin = ^ii'^Qicr'^ajo- (3) 

ija 

H\oc = V{c\,^^C2^a + 4»^Ci,^) -I- JSi,S2i 
2 

+ X! [(^" ~ M)cL<TCm(T + UaUai^naii (4) 

where a = 1, 2 is the orbital index, Sai the spin at site i 
in band a and Jicio- the number of spin a electrons at site 
i in band a. The model ^ and various limits of it (which 
include the Anderson lattice or Kondo lattice, models for 
bilayer compounds or multi-orbital systems) have been 
widely studied. 

The main motivation why we are interested in the 
model ^ in the context of this paper (rather than e.g. 
the one-band Hubbard model) is that it has a simple, al- 
most trivial limit which allows us to discuss the validity of 
Luttinger's theorem for Mott insulators and crossovers be- 
tween different types of insulators. We consider the limit 
where (i) the system is half-filled as there are two electrons 
per unit cell, (ii) the ground-state of Hl^^ is a singlet sepa- 
rated by a gap A from all excited states (see Appendix A), 
and (iii) the hopping is very small, t <^ A. The physics 
of this limit can be fully understood using a straightfor- 
ward strong coupling expansion around the atomic limit, 
Hkin = (performed e.g. up to 11th order in Ref. ^ for a 
Kondo insulator). In contrast, the Mott insulating phase 
of the single-band Hubbard models is more difficult to an- 
alyze as the ground-state of -ffioc is spin-degenerate which 
usually leads to some form of magnetism (and also the 
zero-temperature non-magnetic state obtained e.g. within 
dynamical mean field theory is non-generic and responds 
in a singular way to tiny magnetic fields 7J). Within our 
model one can easily study the crossover from ei = e2, 
V — Q and Ua ^ where one has two coupled Mott 
insulators with A — 3J/4 and charge gap U/2 + 3J/4 
to a band insulator with C/q, = and a large band gap 

A = \J (^^^)^ + V^. Taking a slightly different limiting 
procedure [U2 = 2{n - £2) 00, if^ ^ Q^Ui = 0,V = 
and J > 0], the model ([2]) also describes a Kondo insula- 
tor, i.e. the insulating phase of a half-filled Kondo-lattice 
model. 



2 Breakdown of Luttinger's theorem 
2.1 Local limit 

First, two coupled Mott insulators are analyzed with V ~ 
0, Ui = U2 = U, 61=62 = -u/2, = tl = Uj, J, U » 
tij and /I chosen such that the system is half filled, n = 2. 



For tij = 0, the energy (per site) of the singlet ground- 
state is — 2/i — 3 J/4. The T = Q Greens function of each 
of the two bands is given by 



Gioc(c^) = - ( — = \ — = — 

2\uj + fi-U/2 uj + ii + U/2 



(5) 



with il = U + 3 J/2. The corresponding local self-energy 
is given by 



^loc(t^) 



ii {U/2f 



/I 



(6) 



This form of the self-energy is well known from the Hub- 
bard I approximation [lOj which, indeed, correctly de- 
scribes the zeroth order strong-coupling expansion (note, 
however, that we are not using the Hubbard I approx- 
imation in the following but instead a controlled strong- 
coupling expansion). To leading order in tij, the self-energy 
stays local and the Greens function is obtained as 



G(a;,p) 



(7) 



where the single-particle dispersion ip is the Fourier trans- 
form of tij . To this order, the spectral function at momen- 
tum p consists of two delta-peaks approximately located 
at ±f -I- ^ - /I with weight \ 

It is easy to see that Luttinger's theorem is not valid 
for the Mott insulator discussed above. When the chemi- 
cal potential is changed within the gap, E- < fj, < , the 
total density of electrons, n = 2, on the left-hand side of 
Eq. ll]) remains constant at T = 0. Therefore the choice of 
the chemical potential is completely arbitrary in a grand- 
canonical ensemble at T = 0, (see Sec. 12.31 for a discus- 
sion of a canonical ensemble with fixed particle number). 
At the same time, the right-hand side of Luttinger's sum 
rule, Eq. ([T]), changes from 4 to and is not constant, 
see Fig. [TJ Using for example the leading order strong- 
couphng expansion, Eq. ([7]), one gets G[ijJ = 0,p) > 
for all momenta if < /i < and G{bj = 0, p) < for 
< < E+ where 



E^ 
E_ 



_U_ maxk[tk] 
2 2 

U miuk Uit] 
1 + 

2 2 



O 



O 



J' u 



(8) 



are the lower and upper edges of the spectral gap obtained 
from Eq. ([7]) [due to the finite gap in the system, the strong 
coupling expansion is well behaved, see next subsection]. 

It is well known that for the Hubbard-I approxima- 
tion, the Fermi volume of the doped Hubbard model is 
not constant in violation of Luttinger's theorem (see e.g. 
Ref. [H]). This is an artifact of the Hubbard-I approxima- 
tion. In contrast, we perform a controlled strong-coupling 
expansion, which shows (together with simple general ar- 
guments, see also Sec. |4]) that, the volume within the Lut- 
tinger surface of the undoped system is not fixed by the 
particle number. 
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Fig. 1. Scliematic plot of the Luttinger volume [right-hand side 
of Eq. IT}] as a function of the chemical potential ^. As ^ is 
varied between the lower and upper edge of the gap, E- < < 
E+, the density of electrons, n = 2, remains fixed. While for a 
band insulator (dotted line) Luttinger's theorem (n=Luttinger 
volume) is valid, this is not the case for strong interactions, 
where it is either fulfilled only for a single value of /i, jj, = hl, 
(dashed line) or only a limited range of chemical potentials 
(solid line). For fixed n and T — > 0, the chemical potential 
takes the value fin = {E- + E+)/2 which coincides with fiL 
only for a particle-hole symmetric system. 



are approximately located at 

(9) 

where is the vector pointing from site to site i. By 
definition, the Luttinger surface is given hy fi = E* and 
its volume depends obviously strongly on the value of the 
chemical potential /i in a regime where n = 2 = const, as 
is shown in Fig. [T] As discussed above, Luttinger's sum 
rule is therefore violated. 

Eq. Q is consistent with results of Pairault, Senechal 
and Trcmblay [12] who performed a strong-coupling ex- 
pansion for the single band Hubbard model [in combina- 
tion with a resumniation based on a continued-fraction 
expansion of G(p, uj)] . In their expression, the J of Eq. ([9]) 
is replaced by the temperature T reflecting the degener- 
acy of the strong-coupling ground-state of the single-band 
Hubbard model. 



2.3 Canonical ensemble and position of the chemical 
potential 



This proves that Luttinger's theorem is not valid for a 
Mott insulator. One may, however, ask the following ques- 
tions which will be addressed below: (i) What happens if 
higher orders in the strong-coupling expansion are con- 
sidered such that the self-energy acquires a momentum 
dependence? (ii) What happens in a canonical ensemble 
where the chemical potential is fixed by taking the limit 
T ^ at fixed particle number n? (iii) Which assumption 
underlying the proof of Luttinger's theorem is not valid? 



2.2 Finite bandwidth 

It is convenient to discuss an expansion in the hopping tij 
using the Greens function in real space, Gij^uj), where i 
and j are site indices. To order (ty)^ only the local Greens 
function, Ga gets non-trivial corrections and therefore the 
self-energy remains local. The first correction to Gij 
which is not captured by d?]), arises to order {Uj Y when 
an electron correlates the sites i and j by three hopping 
processes. Close to the poles of G'ioc(w), at ~ ±1] 12, 
one has to resum perturbation theory in to infinite 
order [to obtain e.g. the next correction to Eq. [8])] but 
this is not necessary in the center of the gap where the 
zeros of the Greens function are located. 

Simple power-counting is not sufficient to get the pre- 
cise dependence of the (tijf' contribution on U and J. 
In the limit ^ J ^ C/, by solving a two-site system 
exactly we obtain Gij(w + /i = 0) w — 12(tij)^/(J7^ J). 

Comparing this to Gioc('^) ~ (t^^ + /^)/(C//2)^ for small 
w + /i, one finds after a Fourier transformation that the 
zeros of G(p,w) [and therefore also the poles of i7(p,w)] 



Within the Mott gap, there is precisely one value of the 
chemical potential fi — fj.^ where Luttinger's sum rule, 
Eq. ^ is fulfilled (using again the parameters of Sec. l2.1l in 
this section). Furthermore, in a canonical ensemble (i.e. for 
fixed particle density n) the chemical potential has a well 
defined limiting value for T — > 0, /i„ = liiaT^o ^J-{n,T). 
Therefore the question arises, whether the Luttinger the- 
orem is valid for a Mott insulator if a canonical ensem- 
ble is used where /i = /z„ at T = 0. Note that in a 
three-dimensional system, the electron density is fixed by 
the long-ranged Coulomb interaction and the background 
charge of the ions, therefore the chemical potential does 
take the value /i„ for T ^ 0. Before showing that the 
Luttinger theorem is not valid in this case as generically 
fJ'L 7^ fJ-m one should first note that in the classical proofs 
of Luttinger's theorem [8l[T] always grand-canonical en- 
sembles (fixed /i) and never canonical ensembles (fixed n) 
are used and therefore there is little reason to believe that 
the Luttinger theorem is only valid for a canonical ensem- 
ble. 

For a particle-hole symmetric situation (often studied 
in literature [2,3,^,'^), /i„ and /i^ coincide by symmetry: 
the symmetry fixes both the value of /i and the position 
and shape of the Luttinger surface completely. Generi- 
cally, particle- hole symmetry is not present. To show that 
fj,n and /iL differ from each other in this case, we calcu- 
late both of them to leading order in the strong-coupling 
expansion. 

In the canonical ensemble, the chemical potential in 
located precisely in the middle of the upper and lower 
edge of the gap, as the number of particle excitations oc 
Q-{E+-iJ.)/T Y^g^ equal to the number of holes cx 

g-(M-S-)/T y (see also Fig.[l|). Using Eq. ^ one 
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therefore obtains 



-E+ + max[tp] + min[ip] 



O 



(10) 



This has to be compared to Eq. ([9|) which shows that fiL = 
0+O[{tij)^ / (UJ)]. We have therefore proven that 7^ A'n 
in the absence of particle-hole symmetry when max[ip] 7^ 
— min[ip]. Accordingly, the Luttinger theorem is violated 
in the Mott insulating phase even for a canonical ensemble 
(in the absence of particle- hole symmetry). 

Note that ambiguities in the definition of the chemical 
potential and in the frequency, where the Greens function 
matrix is evaluated within Luttinger's theorem, are ab- 
sent in systems where a Fermi surface exists. Adding to 
our insulating model ^ a further metallic band, which 
hybridizes only weakly with the two bands of the Mott 
insulating state, does not change any of our conclusions 
but fixes /i unambiguously even at T = 



2.4 Analysis of the proof 

We will now investigate the question, which assumption 
underlying the proof of Luttinger's theorem is not valid 
for a Mott insulator. The central step of the proof 0[T] 
(sketched in Appendix B) is to show that the integral 



(11) 



vanishes. The basic idea is to use that the self-energy can 
be written as a derivative of the Luttinger- Ward functional 
with respect to the Greens function such that the relevant 
integrand can be identified with a total derivative. Note 
that in (fTTj) the integrand drops rapidly with l/cj^ for 
large frequencies and that (with the possible exception 
of the point uj — 0) no singularities are expected on the 
integration path along the imaginary axis. As Luttinger's 
theorem is violated in the local limit, one can use the 
exact Greens function and self energy for = 0, given 
by Eqs. ^ and ([6]), to calculate (e.g. by a trivial contour 
integration) 



^ Z'loc(«t^)-^Gloc(«t^) 



27r 



dlLU 



1 

-sign^ 



(12) 



This well-behaved integral does not vanish, contrary to 
what is suggested by the proof of Luttinger's theorem. 
Taking into account spin- and orbital summations, the 
factor ±1/2 explains the violation of Luttinger's theorem 
by ±2 for 7^ in the local limit. 

It is presently not clear, where precisely the problem 
is located in the arguments which try to show that (|lip 
vanishes (see Appendix B) . Despite the simple form of the 
self energy ([S]), bare perturbation theory is very singular 
for /i 7^ 0, r = 0, as the density of the half-filled model, 
n = 2, is reached only above a critical value of U. Inspec- 
tion of the exact finite temperature Greens function in the 
local limit indeed shows divergencies with powers of l/T 




Fig. 2. Evolution of the energy levels in the local limit, tfj = 0, 
see Appendix A, as a function of a control parameter A. For 
A = 0, the inter-orbital hybridization V vanishes, V = 0, and 
Ui = 1,U2 = 1.5, p = 0, J = 0.2, ei = -0.6,62 = -0.5. For 
0<A<1,V = A increases linearly with A and the in- 
teractions and fj. stay constant. For l<A<1.5,V = lis 
constant [15] but the interactions Ui, U2 and J drop propor- 
tional to 4(1.5 — A)^. For A = 1.5, the non-interacting sys- 
tem is characterized by two single-particle levels with energies 
£i±£2 =p y^(£i^£l)%ry2 ~ _i.55 and 0.45. The groundstate 
is always a singlet with two electrons per site, n = 2 well sepa- 
rated by a finite gap from all excited state. As there is no level 
crossing, the wave function evolves smoothly from a singlet 
made from localized electrons in two different bands (Mott in- 
sulator for small A) to two electrons filling a single band (band 
insulator at A = 1.5). The finite gap to excited states guar- 
antees that also in the presence of small but finite hopping 
tfj two-orbital Mott- and band-insulator are adiabatically con- 
nected along the described path. 



to arbitrary order. One may therefore suspect that the re- 
placement of Matsubara sums by integrals, i.e. the limit 
T — > 0, is problematic. Similarly, it may not be sufficient 
to show that (jlip vanishes order by order in a skeleton 
expansion (see Appendix B) for such a non-perturbative 
problem. It is, however, interesting to remark, that the 
Luttinger- Ward functional can in principle be constructed 
non-perturbatively as pointed out by Potthoff [13] who 
also emphasized in his paper that the proof of Luttinger's 
theorem requires non-trivial assumptions on the regularity 
of the limit T — > 0. A more detailed discussion of assump- 
tions underlying the proof can be found in Appendix B. 

In Ref. [M], Altshuler et al. studied the question, to 
what extent Luttinger's theorem applies to an antiferro- 
magnetic metal if one insists to describe the system with- 
out explicitly breaking the symmetry. The authors found 
that in this situation the integral pT|) does not vanish. In 
their case, the problem could be traced back to an anomaly 
(i.e. a contribution arising after the proper regularization 
of a divergence) which seems not to be the case for a Mott 
insulator with unique ground state studied here. 
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Fig. 3. Both eigenvalues of the Greens-function matrix (solid 
and dashed) as a function of u + fi in the local limit tfj — 
for V = 0.1, Ui = 1,U2 = 1.5, J = 0.2,61 = -0.6, e = -0.5 
and a chemical potential within the gap, iJ_ < fj, < E+. Both 
Greens functions change their sign within the gap for this value 
of V. Luttinger's theorem holds only if the chemical potential 
is located between those zeros, see middle panel of Fig. |4l 

3 Adiabatic continuity 

Based on the results discussed above, we will now inves- 
tigate the crossover from a band- to a two-orbital Mott 
insulator (or Kondo insulator) with two electrons per unit 
cell, n = 2. Note that we are not discussing the stan- 
dard single-orbital Mott insulator with an odd number of 
electrons per unit cell, which in d > 1 usually becomes 
magnetically ordered for T = 0. For a band insulator Lut- 
tinger's theorem is always trivially fulfilled and there is no 
Luttinger surface for any value of fi. In contrast, we have 
shown that for V,tij ^ U, J, Luttinger's theorem is vio- 
lated for almost all values of /i within the gap. How does 
the crossover happen or is there a quantum phase transi- 
tion between these two qualitatively different states? This 
question is answered by increasing yv^ + (ei — £2)^ in 
the model ^ in the trivial limit of small tij . 

We first investigate the evolution of the ground state 
when changing the parameters smoothly starting from two 
coupled Mott insulators {V = 0, Ui = l,U2 = 1.5, J = 
0.2) and ending with a band insulator (V — 1, Ui — 0, J — 
0). For tfj = 0, it is shown in Fig. [2] that on this adiabatic 
path the ground state remains a singlet always separated 
by a finite gap A from all excited states (as guaranteed by 
level repulsion). As the gap is finite, these conclusions hold 
also in the presence of a small but finite band-width, tfj 7^ 
0. This proves that there is no quantum-phase transition 
[Is] , separating the two-orbital Mott insulator from the 
band insulator, instead there is just a smooth crossover 
(consistent with results in the literature, see e.g. Ref. [T6]). 

As an aside, we also note that it is generally believed 
that the magnetically ordered states of the particle-hole 
symmetric, half-filled one-band Hubbard model at small 
U (spin-density wave insulator) and large U (magnetically 
ordered Mott insulator) are adiabatically connected. 




V 

Fig. 4. Crossover from the Mott insulator (small intra-orbital 
hybridization V) to a band insulator (large V). For a Mott in- 
sulator with a chemical potential in the shaded region, Lut- 
tinger's theorem is violated (with the exception of the line 
/X = in the upper panel). Luttinger's theorem is, however, 
valid for a band insulator at large V. Solid lines denote the 
upper and lower band edge, the dashed lines the positions of 
the poles of the self-energy or equivalently of the zeros of the 
two Greens functions (see Fig. [3J obtained from the exact so- 
lution of the local model, 1"^ = 0, sketched in Appendix A. 
In the plot the shaded region is not exactly delimited by the 
dashed lines to take schematically a small but finite hopping 
tfj into account (c.f. Fig. [l|. The dotted lines denote the 
position of the chemical potential for fixed particle number, 
fi„ = liniT^o IJ-{n,T). Upper panel: particle-hole symmetric 
case with Ui = U2 = 1, ei = £2 = -1/2, J = 0.2, tfj = 0. 
Middle panel: same parameters as in Fig. [2] Lower panel: 
Ui = l,U2 = 2.5, ei = -0.5, £2 = -2, J = 0.2, t° = 0. 



To investigate the role of Luttinger's theorem and the 
Luttinger surface as a function of the band splitting given 
by ^/V^ + (e2 — ei)^, we first consider the effect of a po- 
tential Ae ~ €2 — El which shifts the two bands rela- 
tive to each other, ei/2 = —U/2 ± Ae/2, in the limit 
V = tfj = and Ui ^ U2 ^ U. The total charge gap 
is given hy U — Ae and the zeros of the Greens functions 
of the two orbitals split by Ae. For Ae < U/2, the two 
zeros remain within the gap. While Luttinger's theorem 
is valid for —Ae/2 < /i < Ae/2, it is violated outside of 
this regime as is shown in Fig. [T] (for tij = the two 
jumps in the solid line become sharp step functions). For 
U/2 < Ae < U, Luttinger's sum rule ^ is valid for all 
values of /i within the gap but the system is still a Mott in- 
sulator in the sense that each of the two bands is singly oc- 
cupied. For Ae = U, two levels cross and the ground-state 
wave-function changes in a first order transition which is 
an artifact of the limit V = 0. Therefore this limit does 
not describe a smooth transition from a two-orbital Mott- 
to a band insulator, which can, however, be obtained by 
increasing the inter-orbital hybridization V. 



6 



A. Rosch: Breakdown of Luttinger's theorem in two-orbital Mott insulators 



Increasing V has a very similar effect as increasing Ae 
for Ae < U a,s is shown in Figs. [3] and ID The shaded 
regions in Fig. |4] indicate the range of chemical potentials 
where Luttinger's theorem does not hold. As discussed 
above, also a finite but small hopping tfj does not change 
the picture qualitatively (see Figs. [1] and [4]) . 

4 Discussion 

The main message of this paper is negative: Luttinger's 
theorem in the form given by Eq. ([T]) is not valid in a Mott 
insulator and therefore the concept of a Luttinger surface 
(points in momentum space where the Greens function 
vanishes at w = 0) seems not to be very useful — especially 
when compared to the much more powerful concept of a 
Fermi surface. 

The basic argument is that a change of the chemical 
potential within the charge gap, E- < fj, < Ej^, just shifts 
the frequency argument of all Greens functions 

G{Lu,p,^^) = G{u; + ^l,p) (13) 

If the Greens functions change their sign within the gap 
(as is the case for a Mott insulator), this implies that the 
sign of G{uj = 0,p) and therefore the Luttinger volume 
depends on the arbitrary choice of fj, while the density of 
electrons remains constant. We have also shown, that even 
in the canonical ensemble, when the chemical potential 
is fixed to limT-*o 2^) by the particle density n, the 
Luttinger theorem is violated in the absence of particle- 
hole symmetry. 

Often Luttinger's theorem is stated in a different way: 
the volume enclosed by the Fermi surface (defined by poles 
rather than sign changes of the Greens functions) is given 
by the number of electrons per unit cell modulo 2 (to take 
into account filled bands). For this case alternative proofs 
exist [l8lfT9] which are based on Gauge invariance, topo- 
logical arguments and the assumption that a Fermi liquid 
exists |18) or on particle-number conservation and an adi- 
abatic connection of the interacting to the non-interacting 
metal (TH]. As we have consided only insulators with n = 2, 
Luttinger's theorem in this "Fermi-surface version" is triv- 
ially valid. Recently, Senthil, Sachdev and Vojta [TT] ar- 
gued that one can construct a situation where this is not 
the case. They considered a two-band model, where the 
first band is in its Mott insulating state such that it can 
be described by free spins. Senthil, Sachdev and Vojta 
then assumed that one can add strongly frustrating in- 
teractions between those spins such that the spins form a 
gapped spin liquid without breaking any symmetry. If the 
second band is weakly interacting, it will form a Fermi 
surface. As the spin-liquid is gapped, this state is robust 
against small perturbations like a finite coupling to the 
Fermi liquid of the second band, e.g. by a small V. In such 
a situation, the volume enclosed by the Fermi surface is 
given by the number of electrons in the second band n2 
rather than by the total number of electrons per unit cell, 
n = 1 + ri2 . Combining this argument with the results of 
our paper, one realizes that in such a case both versions 



of Luttinger's theorem are not valid (in the absence of 
particle- hole symmetry). 

Luttinger's theorem can also break down, if there is 
a finite imaginary part of the self energy at a; = 0. A 
large imaginary part has for example been observed in nu- 
merous studies of weakly doped Mott insulators at small, 
but finite temperatures, see e.g. Refs. [20], or in almost 
magnetic systems, e.g. Ref. '^T. Also in a slightly doped 
tJ model a violation of the Luttinger theorem has been 
observed in small systems [22]. It may be worthwhile to 
emphasize that the breakdown of Luttinger's theorem dis- 
cussed here is of different origin and probably unrelated. 

Our analysis of Luttinger's theorem, Eq. H]), appears 
to imply a qualitative difference between a band-insulator 
where it is valid and a two-orbital Mott insulator where it 
breaks down. Based on the "phase diagrams" of Fig.[4]one 
could define a "critical" hybridization Vc, e.g. by demand- 
ing that Luttinger's theorem is valid for all V > Vc- We 
have, however, shown that at this "critical" hybridization, 
the ground-state energy (and essentially all observables 
with the exception of the Luttinger surface, the Luttinger 
volume and integral (fTT]l ) is analytic in V — Vc- no phase 
transition takes place. Instead there is a smooth crossover 
from a two-orbital {n — 2) Mott- to a band-insulator. 
Recently, Konik, Rice and Tsvelik [3] investigated cou- 
pled Mott ladders constructing explicitly a state where 
simultaneously a Luttinger surface and several Fermi sur- 
faces (due to small electron and hole pockets) are present. 
While we have not studied such a situation, in analogy to 
the results presented here, we speculate that the system 
studied in Ref. [3] is adiabatically connected to some non- 
interacting model where band-structure effects lead to the 
formation of particle and hole pockets. 

While we have shown in this paper that band- and two- 
orbital Mott insulators can in principle be adiabatically 
connected, one should keep in mind that for many models 
there will be a series of first and/or second-order phase 
transitions to various metallic, insulating and/or symme- 
try broken states when the interactions are increased. The 
best studied case in this context is probably the ionic Hub- 
bard model, see Refs. [23,24,25,26 and references therein. 
Note, however, that the symmetries of this special model 
ensure [23], that the Mott insulating phase has no spin 
gap in contrast to the case investigated in this paper. 

For the future, it would be interesting to nail down 
more precisely under what conditions Luttinger's theorem 
breaks down and which of the assumptions used in var- 
ious proofs is most fragile. This is especially important, 
as Luttinger's theorem can be a powerful tool to classify 
topologically different states of matter at least for metallic 
systems. 
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Appendix A: Spectrum of H^^^ 

The spectrum of the local Hamiltonian is given by the 
following energies: The empty site, n = 0, has vanishing 
energy, Eq = 0, while one obtains E4 ~ 2(ei + £2) + Ui + 
U2 — 4/i for n — A. The energies of the four single-particle 



states, n = 1, are Ei± = {^^) ± \J {^^) + - ^, 
while the 4 three-particle states have Ej,± = ei + 62 + 
[^i+l2MIi±U2^j ± y^( n+c/i-.2-t ^)%y2_ 3^, The energy 

of the three n — 2 triplet states is given by E2t = ei + 
62 + J — 2/i. Furthermore there are three singlet states for 
n = 2 with energies -E2SI < -E'2s2 < £-233 which can be 
determined by diagonalizing the 3x3 matrix 

'2{ei-fi) + Ui V2V 

2(e2 -n) + U2 V2V 

V2V V2V ei +62 -2/i -3 J/4, 

(A.l) 

The analytic formulas for E2si can be obtained by solving 
a cubic equation. As they are not very instructive and 
rather long, they are not displayed here. We only note 
that level repulsion ensures that the three energies are 
not degenerate for V ^ 0. 

The total size of the gap for charge excitations of the 
local model with ground-state n = 2 is determined by 
El- + E3- — 2E2si- For y = it is for example given by 
I J -I- min[ei + f/i, £2 + U2] - maxffi, £2]. 

From the eigenfunctions and eigenvalues of the Hamil- 
tonian, one can easily obtain the Greens function. 



Appendix B: A proof of Luttinger's theorem 

In this appendix we briefly sketch a version of the "proof" 
of Luttinger's theorem. The first step is to insert 1 = 
-^{i^ ~ £k) in the expression for the density of electrons 



using imaginary frequencies [T] at T = 



2T 



(B.l) 



p(2vr 



2e*'^' d 



- which we have shown to be not the case for a generic 
Mott insulator. 

Let us nevertheless try to "prove" it. The self energy 
can be obtained by varying the Luttinger-Ward functional 
<?[G], S{p,iuj) = sJL^) ■ The Luttinger-Ward func- 



tional is given by the sum of all skeleton diagrams. By 
the equation <?[G] = J^^G{p,iuj)'pG{p,i^) we define the 
functional $g with 

Eip, zoj) = <?g(p, iLo) + / G(p', ,,.')^^^^^B.3) 

Consider the following integral over a total derivative 







/ -^[G(p,z^)^>g(p,»c.)] (B.4) 



By inspecting the self-energy style diagrams contributing 
to ^g(Pj«i^) one realizes that a change of the external 
frequency w can be absorbed by a change of the frequency 
of all internal lines and therefore 



-^<pG(P,iw) = / 



(5j>G(p,»^) d 
(5G(p', iuj') diu!' 



;G(p',zu;') (B.5) 



One obtains 



= 



f d 

/ — [G(p,ia;)^>G(p,iu;)] 



G 



diuj 

S<pGiP,ii^) d 
^,p, dG{p',tuj') diuj' 



f d 

/ ^(P,*w)— G(p,2tj) 



(B.6) 



where /^^ — duj J^,t ^''P describes the integra- 
tion of frequencies along the imagninary axis and of mo- 
menta over the first Brillouin zone. For the last equa- 
tion, we have performed a partial integration using that 
lmii^^±ao -S'(p, iw)G(p, iuj) = 0. Using the usual analytic 
continuation arguments [51[T] and the vanishing of ImZ'(p, lu 
0) one can identify the first term in (jB.2p with the right- 
hand side of Eq. ([T]) [to keep the presentation simple, we 
omit the discussion of band indices [51[T]]. Luttinger's the- 
orem is therefore valid if the second term in (|B.2p vanishes 



where in the last step we have exchanged the lu, p with the 
uj', p' integration, renamed the variables and used Eq. (|B.3[) . 
If all arguments leading to (|B.6p arc valid, then Luttinger's 
theorem is proven by Eq. (IB.2p . 

However, in the main text, we have shown by an ex- 
plicit example that (jB.6P does not vanish. Therefore one 
of the steps in (|B.6p is not valid for a Mott insulator. The 
first question is whether all integrals are well defined. This 
seems to be the case as for n there is no singularity on 
the imaginary axis and at least order by order in the skele- 
ton expansion each term vanishes with for w —* ±00. 
Also exchanging the order of integration is probably valid 
order by order in the skeleton expansion, as all integrals 
seem to be absolutely converget. As far as we can see, only 
two problems remain. Either the skeleton expansion used 
frequently above is highly singular or the problem is hid- 
den in the limit T ^ and the use of integrals (instead 
=ef Matsubara sums) on the frequency axis and the use of 
derivatives with respect to frequencies was not correct. In 
the main text, it is argued that bare perturbation theory is 
very singular for T — )■ 0, indicating that the two problems 
are related. 
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